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Abstract 



The scaling argument developed by Larichev and Held (1995) for eddy amplitudes and 
fluxes in a horizontally homogeneous, two-layer model on an /-plane is extended to a /3-plane. In 
terms of the non-dimensional number £ = U/(j3\ 2 ), where A is the deformation radius and U is 
the mean thermal wind, the result for the RMS eddy velocity V, the characteristic wavenumber 
of the energy-containing eddies and of the eddy-driven jets kj, and the magnitude of the eddy 
diffusivity for potential vorticity D, in the limit £ >> 1, are as follows: 

v/u « £ ; kjX « r 1 ; « £ 2 

Numerical simulations provide qualitative support for this scaling, but suggest that it underesti- 
mates the sensitivity of these eddy statistics to the value of £. A generalization that is applicable 
to continuous stratification is suggested which leads to the estimates: 

V « (/3T 2 )" 1 ; % « /3T; D « (/3 2 T 3 ) -1 

where T is a time-scale determined by the environment; in particular, it equals XU^ 1 in the 
two-layer model and N^fdzll) -1 in a continuous flow with uniform shear and stratification. 
This same scaling has also been suggested as relevant to a continuously stratified fluid in the 
opposite limit, £ <C 1. (Held, 1980). Therefore, we suggest that it may be of general relevance 
in planetary atmospheres and in the oceans. 



1 Introduction 



An understanding of the mechanisms which control the amplitudes of baroclinic eddies and 
their transport properties is basic to developing theories for the general circulation of planetary- 
atmospheres and oceans, and for incorporating the effects of mesoscale eddies in large-scale ocean 
models. Different classes of idealized models are needed to improve our understanding of various 
aspects of this problem. We focus on one such class - models of horizontally-homogeneous quasi- 
geostrophic turbulence, with imposed large-scale potential vorticity gradients. These models 
are useful in bridging the gap between studies of homogeneous two-dimensional turbulence and 
studies of baroclinic eddy fluxes in inhomogeneous flows of physical interest. Analyses of these 
baroclinic homogeneous models, and work on closely related problems, can be found in Rhines 
(1977), Salmon (1978, 1980), Haidvogel and Held (1980), Vallis (1983), Hoyer and Sadourny 
(1982), Panetta (1993), and Held and 0'Brien(1992). 

Recently, Larichev and Held (1995) (LH hereafter) have considered the statistically steady 
state of a two-layer model on an /-plane, with an imposed environmental vertical shear, or 
interface slope, that results in equal and opposite potential vorticity gradients in the two layers. 
Following the lead of Rhines, Salmon, and Hoyer and Sadourny, they describe a picture of the 
barotropic and baroclinic energy cascades that leads to simple scaling arguments for the energy 
level and potential vorticity fluxes in this baroclinically unstable system, given the scale to which 
the inverse energy cascade extends. 

On a /3-plane, the barotropic inverse energy cascade is halted when the characteristic 
overturning time for the energy-containing eddies becomes comparable to the inverse of the 
Rossby wave frequency (Rhines, 1975), at which scale the energy is channeled into zonal jets. 
In a baroclinically unstable flow, these jets also organize the instability, so that the eddies form 
a 'storm track' on each jet, as described by Williams (1979) and Panetta (1993). By combining 
the scaling argument of LH with the Rhines scale at which the inverse cascade halts, we obtain 
a qualitative theory for the eddy amplitudes and fluxes on a /3-plane. We discuss the results 
for the two-layer model, describe some numerical solutions that provide partial support for the 
theory but also point to some deficiencies, and then indicate how the argument generalizes to 
continuous stratification and arbitrary vertical structure in the mean flow. 

2 The two- layer model 

We consider a two-layer QG model with equal-depth layers when at rest. The barotropic and 
baroclinic streamfunctions are denoted by ip = (^i + ^2)/ 2 and r = (^i — V ; 2)/2, where the 
subscripts 1 and 2 refer to the upper and lower layers respectively. The corresponding barotropic 
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and baroclinic velocities are referred to as (u^, v$) and (u T , v T ) The mean flows in the two-layers 
Ui = const are specified to provide a vertical shear 

U = (U!- U 2 )/2 > (1) 

which results in the mean potential vorticity gradients (5 ± U\~ 2 = £/A~ 2 (l ± where A is 
the internal radius of deformation and where 

e = u/(p\ 2 ), (2) 

The plus sign refers to the upper layer and the minus sign to the lower layer, The lower layer PV 
gradient is reversed if the shear is large enough that £ > 1, which is the criterion for instability 
in an inviscid flow. 

The argument in LH for the /-plane can be reformulated as follows. Start by assuming 
that the kinetic energy of the flow is predominately barotropic on scales much larger than the 
internal radius of deformation A. On these scales the barotropic mode evolves as in a two- 
dimensional flow, cascading energy to larger scales. Assume that this inverse energy cascade 
halts on average at the wavenumber ko- Given the Kolmogorov energy spectrum of A; -5 / 3 , or 
any spectrum steeper than A; -1 , the bulk of the energy will be contained on the scale k . The 
baroclinic potential vorticity is advected by this nearly barotropic flow on large scales. Since it 
does not induce a significant part of the flow by which it is advected, the baroclinic potential 
vorticity will behave as a passive tracer, and will be mixed downgradient by the 'turbulent 
diffusion' engendered by the barotropic flow. This mixing will be dominated by the largest 
scales in the flow. Since the magnitude of the ambient potential vorticity gradient is U\~ 2 , the 
typical size of the eddy baroclinic potential vorticity will be 

q'^k^{U\- 2 ) (3) 

so that 

Vq 1 a* -D(U\~ 2 ) (4) 

with 

D sd Vko 1 (5) 

where V is the rms barotropic velocity. Thus, (3) is equivalent to the statement that the baroclinic 
velocities on the energy containing scale (v' T ~ k X 2 q') are of the order of U. 

The baroclinic potential vorticity, acting as a passive scalar, will cascade to smaller 
scales. On scales larger than A, the variance of this potential vorticity is dominated by thickness 
fluctuations, and is therefore proportional to the potential energy in the baroclinic mode. The 
rate of eddy energy production per unit mass in a horizontally homogeneous two-layer model is 

e = -{u^rg + u 2 ^g)/2 = U^/X 2 » Vk^{U/\f (6) 
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In this homogeneous model, the potential vorticity fluxes in the two layers are equal and opposite 
and proportional to the thickness flux. 

The baroclinic eddy energy is produced on the largest scale A;^ 1 since it is proportional 
to the potential vorticity flux. It then cascades down to the deformation radius. Energy cannot 
cascade further since the layers decouple at smaller scales, leaving 2D flow in which energy cannot 
proceed downscale; instead it is converted to barotropic energy, in which form it cascades back 
upscale, eventually to be dissipated. 

The barotropic energy level is determined by the requirement that at equilibrium the rate 
at which energy flows to larger scale in the barotropic mode is equal to the baroclinic energy 
production. If the energy cascade extends to the scale ko and the RMS barotropic velocity is V, 
then by dimensional analysis, 

e « V 3 k (7) 

an expression familiar from 3D turbulence (e.g., Tennekes (1972)). It is equivalent to the as- 
sumption utilized by LH that an inertial range exists. From (6) and (7), one finds 

V « (hX^U (8) 

The barotropic velocities are larger than the mean shear by the factor (/c A) _1 . We have as- 
sumed that this factor is much greater than unity so as to have a significant inverse energy 
cascade. Barotropic velocities are also larger than baroclinic velocities on the scale k by the 
same factor, so the dominance of advection by the barotropic flow on large scales is assured. 
The corresponding diffusivity is of the magnitude 

D « (koX^Uko 1 (9) 

We also have the result that the eddy available potential energy is proportional to the barotropic 
kinetic energy. 

A^(W^ + M " (10) 
The argument above is equivalent to that provided by LH. 

In the basic /-plane numerical experiment analyzed in LH, k$ 1 is determined by the size 
of the domain. On a /3-plane, we expect the inverse cascade to stop at the Rhines scale, at which 
rms barotropic velocities V are comparable to Rossby wave phase speeds 

K 2 = v/p (n) 

We ignore the anisotropy of the dispersion relation in this scaling. By combining (7) and (11), 
one obtains the expression for k in terms of energy cascade rate discussed by Pelinovskiy (1978) 
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or Vallis and Maltrud (1991): ko (/^/e) 1 / 5 . (If the lower boundary is not flat, one should 
replace (3 by the effective f3 felt by the barotropic mode in the presence of large-scale bottom 
relief, H\V(f /H)\). As the cascade is halted, the barotropic energy is organized into jets of 
this meridional scale. Panetta (1993) has shown that this relation accounts quite well for the 
number of jets produced in a model of homogeneous geostrophic turbulence in the presence of 
a mean shear, if V is set equal to the square root of the total kinetic energy (see his Figure 4). 

Our assumption that (/coA) -1 ^> 1 translates into the statement that £ 3> 1. We assume 
that the presence of f3 does not modify the previous scalings when £ is large. Combining (8) 
and (11), we have 

V/U » £ (12) 

and 

koX^C 1 (13) 

with the consequence that 

D « UXi 2 (14) 

Since D is proportional to U 3 , the potential vorticity flux is proportional to U 4 while the energy 
generation is proportional to U 5 . 

3 Some numerical results 

These scaling relations can be tested against numerical experiments such as those described by 
Haidvogel and Held (1980) and Panetta (1992), for 2-layer, QG homogeneous turbulence on a 
/3-plane with imposed mean vertical shear. Figure 18 in Panetta (1993) shows that the scale 
with the maximum energy and the scale of maximum energy generation are both approximately 
linear in £, consistent with (13). (In the notation of Panetta, £ = (2/5) The predicted 
eddy amplitudes and fluxes are not that well verified, however. The eddy energy should be 
proportional to £ 2 . The results tabulated by Panetta suggest that the energies and fluxes are 
even more sensitive to £. To confirm this result, we have obtained numerical solutions to the 
same model as considered by Haidvogel and Held (1980) and Panetta(1993), but with 256x256 
resolution that allows more scale separation between the radius of deformation and the energy 
containing eddy scale. 

The numerical model is identical to that used in the /-plane simulations of LH. We choose 
a value of the lower layer Ekman damping, k = 0.16(C//A) that is identical to that in LH, and 
set the deformation radius so that 27rA = L/50, where L is the size of the square doubly-periodic 
domain. Small scale damping is modeled with a iA7 8 diffusion operator, with v = 0.08(£/A 7 ). 
Statistically steady states have been obtained for the values of £ shown in Table 1. Instability 
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Table 1: Statistics from the numerical model as a function of £: erfj SS is the dissipation of energy 
in the barotropic mode by Ekman friction; k is computed from (11) using the rms meridional 
velocity of the barotropic mode for V. 

actually persists for values of £ slightly less than 1, due to dissipative destabilization, but these 
are not displayed in the table. 

Figure 1 shows the horizontal spectra of the barotropic component of the meridional 
velocity and of the baroclinic eddy energy generation (or equivalently, of the potential vorticity 
flux in either layer, or of the buoyancy flux) for each of the statistically steady states obtained. 
The spectra are plotted as a function of total horizontal wavenumber, after averaging over angle 
in k-space. The energy generation clearly moves to larger scales when the energy is allowed to 
cascade to larger scales. Also shown by arrows is the Rhines scale (11), computed using the rms 
barotropic meridional velocity for V, for each value of £. This prediction of the scale is seen to 
be in reasonable agreement with the scale of the maxima in energy and energy generation, for 
large values of £ down to £ 2. 

As Fig. 2 shows, the energies and the energy generation increase more rapidly with in- 
creasing £ than predicted by our scaling arguments. Instead of the expected £ 2 power law, we 
obtain a dependence that is between £ 3 and £ 4 . This result confirms the impression obtained 
from Panetta's results. We cannot rule out the possibility that even the largest values of £ are 
not yet large enough to attain the asymptotic regime predicted, but this appears to be unlikely 
since there is no trend in the results that suggests that the £ 2 behavior is being approached. 
(The comparison with theory is also made more complicated by the fact that the size of the 
domain is beginning to make itself felt at the largest value of £ used.) This discrepancy is also 
seen in the /-plane results in LH, in which the scale of the energy generation was arbitrarily 
constrained. 

To examine the sources of this discrepancy, we first plot in figure 3 the energy generation 
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e normalized by 

^ 2 W) 1/2 (^) 1/2 (15) 

This quantity equals unity if the meridional velocity and thickness are perfectly correlated. The 
values obtained are close to 0.2 and hardly change except for the largest value of £ examined. 
Therefore, the neglect of changing correlations in the estimate (6) does not appear to be a 
significant source of error. 

Combining (6) with the estimate (11) for the mixing length, we have 

(16) 

Alternatively, combining (7) with (11) yields 

~ (-) 3 (A; A) « (-) 5 /V-V2 = (1 7) 

[/ 3 A -i v 1 ; K u> q u^x- 1 y ' 

These two expressions for e can be combined to yield (12)-(14). We plot eje\ and e/e2 in 
Figure 4a as a function of £. Once again V has been set equal to the rms barotropic meridional 
velocity. Considering that e varies by more than 3 orders of magnitude (see table), the constancy 
of these ratios is encouraging. However, it is apparent that the RHS of (16) underestimates the 
dependence of the energy production on £, indicating that the effective mixing length increases 
somewhat more rapidly with £ than does the Rhines scale. In contrast, (17) overestimates the 
dependence of e on £, indicating that the inverse energy cascade is not as efficient as indicated 
by (7) if one uses the Rhines scale for k . 

In Fig. 4a we have also compared the estimate (17) to the simulated energy dissipation 
by Ekman damping in the barotropic mode rather than to the eddy energy production. These 
are implicitly assumed to be equal, or at least proportional, in our scaling arguments, but they 
differ in the numerical model because of the energy lost due to the subgrid diffusivity and the 
Ekman damping of baroclinic kinetic energy (see Table 1). The dissipation in the barotropic 
mode is a better estimate of the rate at which energy is cascading to larger scales, which is the 
more relevant quantity for comparison with (17). The curve becomes a bit flatter at smaller 
values of £, but most of the discrepancy remains. Although we have compromised our resolution 
at small scales in order to provide room for a substantial inverse cascade, resulting in non-trivial 
energy losses due to the subgrid scale diffusivity, we do not think that a model with higher 
resolution and weaker diffusion would produce qualitatively different results. 

If we use the total (zonal plus meridional) rms eddy velocity for V in (7), while retaining 
the use of the meridional velocity in (11), the plot analogous to that in Figure 4a (dashed line) 
shows a much improved fit, as shown in Figure 4b. This choice is admittedly somewhat arbitrary 
in the absence of a theory for the anisotropy of the energy spectrum, but it does indicate that 
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the anisotropy of the eddies (see the ratio of rms eddy and Vf in the table), which we have 
ignored in our simple scaling, likely plays a role in the discrepancy in Fig. 4a. 

Because (16) underestimates the energy production and (17) overestimates it, these errors 
add to produce the larger errors in the dependence of the energy level and fluxes on £. For 
example, if one arbitrarily replaces (16) and (17) with 



)3/2+a£l/2 (lg) 



then one finds 

(y/U) 2 ^i 5 -5^2/{l-a- 1 ) (20) 

Choosing a = 7 = 0.20 or 0.25, as Fig. 4a roughly suggests, one obtains 5 = 10/3 or 4, consistent 
with the variation of energy with £ in Figure 2. Errors that appear to be small in the individual 
approximations combine to create a large discrepancy in the final dependence of energy on £. 
The underlying reason for this sensitivity is the strong positive feedback inherent in the system: 
as the energy level increases, the length scale of the energy containing eddies increases, which 
increases the mixing length and V, increasing the energy generation. 

Although our qualitative arguments clearly require some modification in order to fit 
numerical results, we feel that they are a useful starting point, and therefore have examined 
how they might be generalized to horizontally homogeneous flows with arbitrary stratification 
and mean vertical shears. 



4 Continuous stratification 

Consider a Boussinesq QG flow with stratification N 2 (z) and mean flow U(z), confined between 
flat horizontal boundaries at z — 0, H. The energy production per unit mass e can be written 
in the form 

_! r H fv'b'dU 



_ e = H L V** < 21 > 

Here v'U is the meridional buoyancy flux, related to the potential vorticity flux by 

777= fd z (N~ 2 VF) (22) 

Consistent with our assumption of horizontal homogeneity, the relative vorticity flux, which 
equals the convergence of the eddy momentum fluxes, is ignored. The mean wind shear is 
related to the mean buoyancy gradient by 

fd z U = -d y B (23) 
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We assume again that the inverse cascade is substantial enough that the flow advecting 
the potential vorticity can be taken as barotropic. The downgradient fluxes can then be set equal 
to a diffusivity D times the negative of the mean gradient, with the diffusivity independent of z. 
In this case of a vertically uniform diffusivity, the same diffusivity can be applied to buoyancy 
and to potential vorticity. The result is 

e = DT~ 2 (24) 

where 

For the two layer model of Section 2, T = X/U, where U = (Ui — U 2 )/2, from which we can 
rederive the 2-layer scaling. For a continuous flow with iV and d z ll constants, we have T = 
N(fd z U)~ 1 = f~ l Ri l l 2 . In both cases, T is proportional to the e-folding time for baroclinic 
waves in the limit that the effect of (3 is small. We again assume that the diffusivity D is related 
to the energy containing scale k and the barotropic velocity scale V, as in (5), and that the 
energy production must balance the upscale transport of energy by the barotropic flow (7), so 
that 

e = Vk^T' 2 = V 3 k , (26) 

or simply, 

T = l/(k V) (27) 

Thus, in addition to being the time scale associated with the linear dynamics of deformation 
scale eddies, T is also the time scale of the energy containing eddies on much larger scales. 
Combined with the estimate of k as the Rhines scale (11), we have 

k = PT, (28) 

V = l/([3T 2 ), (29) 
D = l/(/3 2 T 3 ) (30) 



5 Relationship to more weakly unstable flows 

The preceding analysis is intended for the asymptotic limit of small f3, in which there is an inverse 
energy cascade over a substantial range of scales. It is in this limit that the flow responsible for 
diffusing the potential vorticity is primarily barotropic, so that we can remove the diffusivity 
D from the integral in the expression for the energy generation (21). In the two-layer model, 
these scaling relations cannot remain qualitatively valid for all £, since the flow is stable once £ 
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drops below a critical value (although it is impressive how well the two-layer results in Figure 2 
are approximated by a power law even as £ approaches unity). However, these expressions may 
retain their relevance over the full range of £ in certain continuous flows. Held (1978) suggests 
scaling arguments for Charney-like continuous models in the limit that h/H <C 1, where H is 
the depth of the fluid, or the scale height in the non-Boussinesq case, while 

f 2 d U 

In this limit the vertical scale of the most unstable waves in Charney's model is proportional to 
h, and their horizontal scale to 

Nh/f = ^ = (32) 

where we have used the expression (25) for T. This is identical to the scaling in (28). If 
these scales dominate the statistically steady state, surface buoyancy perturbations will be of 
magnitude b' pa (Nh/ f)d y B, or, using equipartition of eddy kinetic and available potential 
energies, 

V ^ b'/N ^ hd z U pa (3~ l (fd z U / N) 2 (33) 

which is identical to (29). The result is again a diffusivity proportional to the cube of the vertical 
shear. 

In Held (1978), this scaling was assumed to break down as h/H becomes comparable to 
unity. Following Stone (1972), it was assumed that the dominant horizontal scale relevant for 
the eddy transports would then be NH/ f — A and the characteristic velocity Hd z ll, leading to 
a diffusivity that is proportional to the first power of the vertical shear. The argument of LH 
suggests instead that once h/H(= £) rises above unity, the inverse energy cascade sets in and 
the appropriate eddy scale and mixing length continue to increase in proportion to £. With this 
picture in mind, one can imagine that the scalings (28-30) are relevant for the full range of £. 

In the weakly unstable case one cannot assume that the diffusivity is independent of 
height. In fact, if h/H <C 1 the mixing should be confined to a depth proportional to h (Held, 
1978) in a Charney-like environment. More generally, a theory for the vertical structure of the 
diffusivity is required as one moves away from the strongly unstable limit. 



6 Final remarks 

The qualitative theory presented for the eddy amplitudes and scales in horizontally homoge- 
neous baroclinically unstable flows can be most simply described by the statement that the 
environmental shear and stratification determine a time scale T, which can be combined with j3 
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to determine a length (/ST) -1 , a velocity (/9T 2 ) -1 , and a diffusivity (/3 2 T 3 ) -1 in only one way. 
The remainder of the argument provides the rationale for thinking that the shear and stratifi- 
cation provides a time scale only, in the limit in which the flow is energetic enough that there 
exists a substantial inverse energy cascade. The result coincides in form to that proposed by 
Held (1978) for a Charney-like environment, in the opposite limit in which the eddies are very 
weak. This leads us to suggest that these scalings may in fact be useful for a wide range of 
energy levels in certain flows, although we do not understand if there are fundamental reasons 
why this should be so. 

In practice, one must keep in mind that there may be insufficient room for (3 to halt 
the cascade, in which case the eddy scale is determined by the domain size, or the size of the 
baroclinic region. This leads to the diffusivity proposed by Green (1970), as described by LH. 

Several limitations to these arguments are evident. Even in the limit of very strong su- 
percriticality in the two layer model, there is only partial agreement with numerical simulations. 
Our scaling arguments suggest that an accurate theory for eddy amplitudes and fluxes will be 
difficult to obtain, because of the delicate balance resulting from a feedback in which potential 
vorticity fluxes increase as eddy velocity and length scales increase, leading to higher energy 
levels and a stronger inverse energy cascade, leading, in turn, to an increase in eddy velocity 
and length scales. 

Even if the limitations of the asymptotic theory are ignored, one still needs a theory for 
the vertical structure of the diffusivity, or of the eddy PV fluxes, when the inverse energy cascade 
is not so strong as to generate a nearly barotropic eddy velocity field. These arguments also do 
not tell us how to treat the dependence of eddy statistics on the strength of the surface drag. 
Also, the bottom topography so far ignored in the analysis may strongly effect the cascades, 
rendering the flows more baroclinic (Cox, 1985; Treguier and Hua, 1988). 

The final and most fundamental limitation relates to the assumption of horizontal ho- 
mogeneity. For example, if the barotropic decay is due to momentum fluxes in the presence 
of a mean horizontal shear, and if this horizontal shear is imposed by other constraints, rather 
than being internally generated by the eddies themselves, then the dependence of the length 
and velocity scales on environmental parameters would clearly be changed. However, the recent 
study of Pavan and Held (1995) suggests that diffusivities obtained from a homogeneous model 
can be surprisingly useful for jet-like flows with substantial horizontal shears. 

Despite these limitations and misgivings, we believe that these scaling arguments provide 
a useful background for developing closure schemes for simple diffusive atmospheric models as 
well as for non-eddy resolving ocean models. 
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Figure Captions 



Fig. 1. The horizontal spectra of (a) the barotropic meridional velocity and (b) the baroclinic 
eddy energy production rate, for different £ at equilibrium. The arrows indicate the Rhines 
scale computed from (11) using the rms barotropic meridional velocity for V. 

Fig. 2. A log- log plot of v'^ 2 , twice the meridional barotropic energy (solid line), and e, the 
baroclinic eddy energy production rate (dashed line), as functions of the supercriticality 
£. Also shown are slopes consistent with £ 2 and £ 4 dependencies. 

Fig. 3. The baroclinic eddy energy production rate, e as a function of supercriticality £, 
normalized by (15). 

Fig. 4. Ratios of the theoretical estimates (16), (17) to the numerical simulations, as a function 
of the supercriticality £: (a) eje\ (solid line), e/e2 (dashed line), and e<iissl e 2 (dotted line). 
See (16) and (17) for definitions of t\ and t^- tdiss is the energy dissipation by Ekman 
damping in the barotropic mode; (b) e/t2 modified by using the total rms eddy barotropic 
velocity in (17) (rather than the meridional velocity only, as for the dashed line in (a)). 
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Figure 1: The horizontal spectra of (a) the barotropic meridional velocity and (b) the baroclinic 
eddy energy production rate, for different £ at equilibrium. The arrows indicate the Rhines scale 
computed from (11) using the rms barotropic meridional velocity for V. 
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Fig. 2 




Figure 2: A log-log plot of v'^ 2 , twice the meridional barotropic energy (solid line), and e, the 
baroclinic eddy energy production rate (dashed line), as functions of the supercriticality £. Also 
shown are slopes consistent with £ 2 and £ 4 dependencies. 
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Figure 3: The baroclinic eddy energy production rate, e as a function of supercriticality £, 
normalized by (15). 
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Fig. 4 
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Figure 4: Ratios of the theoretical estimates (16), (17) to the numerical simulations, as a 
function of the supercriticality £: (a) e/ei (solid line), e/t2 (dashed line), and edi SS /^2 (dotted 
line). See (16) and (17) for definitions of e\ and e 2 . e^ ss is the energy dissipation by Ekman 
damping in the barotropic mode; (b) e/e 2 modified by using the total rms eddy barotropic 
velocity in (17) (rather than the meridional velocity only, as for the dashed line in (a)). 
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